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Abstract: 

We show that experimental data are in very good agreement with predictions from the string 
fragmentation model by Bowler and Morris. We present a physical interpretation and discuss 
the relation to results obtained from Perturbative QCD and Local Hadron Parton Duality 
(LPHD). We also present implications for B7r correlations and the possibility to use these as 
a tag to study CP violation in B decays. 
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1 Introduction 



An understanding of hadronization in b-quark jets is important for several reasons. 

• It can help to understand the confinement mechanism. 

• A good description of jet fragmentation is essential to reconstruct an underlying hard 
parton interaction. Thus e.g. a measurement of the reaction e + e~ — > W + W~ — > 
q 1 q 2 q 3 q 4 at LEP2 is important for determination of the W mass, which can give infor- 
mation about virtual corrections involving the Higgs particle. 

• As proposed in [|TJ particle correlations in b and b jets can be used as a tag which e.g. 
can improve a determination of CP violation in B decay. 

The Lund string fragmentation model (implemented in the Jetset MC ||) and the 
cluster fragmentation model [|J (implemented in the Herwig MC ||) have been frequently 
used by experimental groups to describe the hadronization process in the analysis of their 
data. Previous official versions (up to v7.3) of the Jetset MC have in the default version given 
a too hard B meson spectrum, and in consequence too low total multiplicity in bb events. To 
describe the data, experimentalists have instead often used the Peterson et al. model || for 
the B meson momentum, combined with string fragmentation for the remainder of the jet. 
That model contains one free parameter which can be adjusted to the experimental data. 
The two most essential experimental observables are the average B-meson energy (or Tb) and 
the average total charged multiplicity n. In the Jetset MC it is assumed that after separation 
of the B hadron, the rest of the jet corresponds to a light quark jet with the remaining energy, 
and therefore there is a correlation between xb and n. It is noteworthy that this correlation 
agrees well with data. Bowler and Morris ]7| have proposed a modification of the initial Lund 
Model for the fragmentation of heavy b or c quarks. This model has no free parameter, and 
therefore gives a definite prediction. In the present paper we will demonstrate that the model 

• agrees well with expectations from a string dynamics scenario. 

• agrees qualitatively with expectations from perturbative QCD and local parton hadron 
duality or cluster fragmentation. 

• agrees very well with experimental data. 

Equipped with a description of b fragmentation, we will in this paper also study the Bit 
correlations and estimate the efficiency of the tagging method proposed in ref M for a study 
of CP violation in B decays. 

2 b jet fragmentation 

The Lund string fragmentation model is implemented in Jetset to take care of the frag- 
mentation of qq strings. The main idea in the fragmentation of the original qq pair, ie. 
e + e~ — > Z°,7 — > qq, is that new quark-anti-quark (diquark-anti-diquark) pairs are pro- 
duced in the colour field stretched between the original quark pair. The constituents of 
a new pair obtain their transverse mass through the tunneling process needed to create 
them ||. There is a flavour ordering present throughout the string so that the flavour quan- 
tum number is preserved. Details of how these properties are given to hadrons produced in 
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Figure 1: Quark trajectories in the Lund model. Light quarks move along light cones a), 
while heavy quarks move along hyperbolae b). The finite mass of heavy quarks decreases the 
area bounded by the solid lines resulting in softer B spectrum. 

Jetset will not concern us here. Instead we are going to study the energy and longitudinal 
momentum distribution of the produced hadrons. 

A basic assumption in the Lund model is "left-right symmetry" , which means that the 
fragmentation process should look the same irrespective to which end of the string we start 
from in an iterative process. From this assumption follows that the probability dP to obtain 
a definite final state with hadron momenta (being 1+1 dimensional vectors if we neglect 
the transverse dimensions) with masses is given by || 

dP oc J] N^p^^pj - m 2 t )5(Y,Pk ~ Ptotai) exp(-bA) (1) 

i k 

where A is the area indicated in Fig [l]. It is suggestive that this expression is the product of 
a phase space factor and the exponent of an area, which can be interpreted as (the imaginary 
part of) an action and is similar to a Wilson loop integral. 

The distribution in Eq ([!]) can be generated iteratively starting from one end of the 
quark-ant i- quark system. With the other properties fixed, the probability to give a hadron 
a fraction z of the available energy is given from a splitting function, 

The parameters dj could in principle depend on the flavour of the associated quark-anti- 
quark pair. In Eq (|) the index i corresponds to the flavour of the previously produced qq 
pair and k to the latest produced flavour pair. Since experimental data seems to need the use 
of only two different a- values in connection with quark-ant i- quark and diquark-anti-diquark 
production respectively (the latter leading to baryon-anti-baryon production), two different 
a-values have been used in most applications. 

We also note that in the string model it is assumed that heavy quarks like c or b can only 
be produced perturbatively, either directly from the initial 7 or Z, or in the process g — > qq, 
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and not in the soft string breaking process. Thus a b-quark is always at the end of a string, 
and the argument based on left-right symmetry is really not applicable for the B hadrons. 
The fragmentation will always be left-right symmetric whatever energy is given to the B and 
the B at the ends of the string. A modification of the splitting function in Eq @, based on 
the physical ideas behind the string model, was proposed by Bowler and Morris in 0. If the 
initial qq pair is a heavy cc or bb pair, these quarks do not move along the light cones but 
instead along hyperbolae in the x — t diagram as shown in Fig (TJb. It is then natural to insert 
in Eq ([[]) the area bounded by the quark trajectory as indicated in Fig [lfb 0. The result is 
a softer spectrum for the leading B meson, which is well approximated by the expression 



Here rq is predicted to be equal to 1, but is introduced as a parameter in Jetset. This is 
done to make it possible to control the effect of the modification. 

A relativistic string or a homogeneous (colour-)electric field is invariant under longitu- 
dinal boosts, which in the string model is reflected in a smooth distribution in rapidity for 
the produced hadrons. Although it is not at all a definite consequence, we feel that from the 
physical picture of hadrons formed from the energy stored in a string-like field, it is most 
natural that the distribution of light hadrons stretch up to the rapidity of the leading heavy 
B meson (or D meson in case of a c quark jet). Thus there should neither be a large rapidity 
gap nor should the rapidity distribution of the light hadrons continue beyond the rapidity of 
the heavy meson. This is actually the case for the splitting function in Eq (|3]), while for the 
function in Eq (Q) a large gap is obtained between the leading B-meson and the remainder of 
the jet. The average rapidity difference, Ay, between the B meson and the first rank meson 
in the remainder is 0.88 which should be compared with the rapidity difference 0.91 between 
neighbouring light mesons in the middle of the string. The corresponding numbers for the 
harder distribution is 1.96 and 0.92, respectively. (As representative for a light meson we 
have used p mesons to avoid kinematical effects from the exceptionally low pion mass, cf 
the discussion in ||.) The difference between the two distributions is further illustrated in 
Fig |2|. This figure shows the rapidity difference distributions between differently flavoured 
neighbouring mesons. The large (compared to the light mesons within the string) rapidity 
difference between the leading heavy B and the light p is present for the hard distribution 
while the softer gives similar distributions irrespective to flavour and rank. 

Such a smooth rapidity distribution is also expected in PQCD assuming LPHD or cluster 
fragmentation. This is a consequence of the fact that the initial b-quark radiates gluons with 
(pseudo-)rapidities up to the rapidity of the b-quark, but not beyond this value. The region 
of larger rapidities has been called the dead cone [JHJ . 

To be able to describe LEP e + e~ — > qq data correctly gluons must be included. In 
all simulations below gluons are radiated before string breakup. When we use the original 
symmetric Lund fragmentation formula (Eq @) in the simulations we find that the average 
xe for B-mesons decreases to 0.77. This is well outside the experimental value, 0.702±0.002± 
0.008 (average of the results reported from L3, ALEPH, OPAL and DELPHI experiments at 
LEP in 1994) JLTJ]. Using the modified formula, Eq (|3|), we find that tq = 1.05 ± 0.07, gives 
an xe within the experimental errors. We note that this is in very good agreement with the 
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Figure 2: a) The full line shows the rapidity difference between rank 2 and rank 3 p's, 
and the dotted line shows the rapidity difference between B and rank 2 p's, obtained by 
the modified splitting function, Eq (||). b) Shows the same distributions when the original 
splitting function, Eq (\^), is used. 
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Figure 3: Combined B°B + energy spectrum obtained at the Z pole, using the original (dotted 
line), and the modified (solid line) splitting functions, Eqs ([Jj and (^j, respectively. (Gluon 
radiation is included before fragmentation.) 



theoretically expected value 1. In Fig ^| we show the combined B° and B + energy spectrum 
using Eq (g) and Eq (^), respectively. We see that the spectrum is considerably softened by 
the modified function. 

The acceptable values of r<g shift down to 0.98±0.08 when 30% of the produced mesons 
have orbital angular momentum (tensor mesons) in the simulations. The B meson energy 
spectrum is softened by the (B7r) resonance states, and in consequence the splitting function 
can be harder as compared with simulations without tensor mesons. 

Using the softer b hadron energy distribution (Eq (|3|) with rg = 1) will increase the 
average number of charged particles produced in Jetset. The average difference between 
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±0.002 
±0.008 
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Table 1: Some typical results obtained when the fragmentation models are tuned to retain 
the charged multiplicity at LEP1. The 'Tensor' column shows the fraction of tensor meson 
production in the simulations, and the 'Corr. ' column shows if flavour correlations are 
included. 



the number of charged particles produced in b events compared to uds events, n\> — n u d s , 
becomes 3.06^o;og (2.63l[j;2g f° r 30% tensor mesons). This is in agreement with the exper- 
imental value, 3.14 ± 0.44 [[LI]]. Compared with an analytic MLLA calculation which gives 



5.5 ± 0.8 [0, there is a large discrepancy. Despite this discrepancy we see that a detailed 
MC, with the same fundamental properties as MLLA, agrees with experimental data. The 
original splitting function in Eq (0), however, gives an average difference between the number 
of charged particles in b and uds events of 1.21, which is well below the experimental value. 

Similar calculations have also been performed using other fragmentation models im- 
plemented in the Jetset MC: the Peterson et al. model for heavy quark fragmentation, 
inclusion of tensor meson production, and introduction of transverse momentum correla- 
tions and flavour correlations between neighbouring string breakups presented in |J. In all 
simulations an optimization of the fragmentation parameters, a, b, and e&, have been per- 
formed in order to try to retain several experimental observables (n c h, — n u d s , x~e for B 
mesons, and average thrust) at LEP1 energies. In Table |1] we list results obtained when 
the different fragmentation models are used. The a and b parameters are chosen to give the 
best compromise for a given model, with the condition that the observed charged particle 
multiplicity at LEP1 is retained. 

We see that the Peterson et al. model, and the original splitting function, do not 
reproduce the experimental data. The modified splitting function performs well irrespective 
if tensor mesons are produced in the fragmentation or not, when no flavour correlations are 
included. However, in the presence of flavour correlations there seems to be a need of tensor 
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meson production. We believe the fraction of tensor meson production used in the MC's are 
extreme values of the real fraction at LEP1, and note that the OPAL Collaboration expects 
the fraction to be at least 20% IlLl. Inclusion of flavour correlations in the simulations 



improve the performance of MC's when the Peterson et al. model or the original splitting 
function is used, but the discrepancy with experimental data is still quite large. 

We conclude that the splitting function in Eq (|3|) is both physically motivated and in 
better agreement with experimental data than the other fragmentation models implemented 
in Jetset, although these cannot be fully ruled out on the basis of above considerations. 



3 B7r correlations as a flavour tag and CP violation 

In it was suggested to use Bn correlations as a flavour tag in a study of CP violation in 
B decays. In this section we will study these correlations and estimate the efficiency in CP 
violation studies. We will also study how the relation between efficiency and purity varies 
with the mass of the B7r pair. 

B7r correlations depend on the properties of b-jet fragmentation, and also on the pro- 
duction rate for B resonances decaying into a Bit pair. In the default version of the Jetset 
MC there are no correlations in flavour or transverse momentum between neighbouring qq 
string breakups. In || it is argued that due to the exceptionally small pion mass (the pions 
have also properties corresponding to a Goldstone boson) such correlations are expected to 
some degree. These flavour correlations could influence the B7r correlations, even if their 
effect on x B and the inclusive spectra is negligible. In the following we will also study this 
possibility. 

If a B° (bd) meson is produced in a bb event the associated d is likely to end up in 
a 7r + near in momentum space. This ir + could be part of the remaining jet, either as a 
directly produced pion or the decay product from e.g. a p meson. The B°7r + pair can 
also be decay products from a parent B resonance. From iso-spin invariance the B°7r + 
correlations are related to B ± 7r correlations, which are experimentally studied by the OPAL 
collaboration |14| . 

Let us first study the case where the d partner of a B° is a leading parton in the 
remaining jet. The tagging possibility relies on the fact that the excess positive charge from 
the d is found in the leading end of the jet, usually in the first rank hadron (of the remaining 
jet). If it is found in the first rank hadron the characteristic invariant Bn mass is determined 
by the relation 

M^ ar ^M B + m ± cosh A' (4) 

where m± is the transverse mass of the pion and A' is the rapidity difference between the 
B and the pion or the parent resonance if the pion comes from e.g. a p meson. Since the 
energy spectrum for the B mesons is strongly peaked at high energies, with an average around 
xe ~ 0.7, the Bit mass increase rapidly for higher rank pions with lower energy in the cms. 
Thus to enhance the signal it is favourable to study only pairs with mass below 5.8GeV. 
We also note that low energy B mesons can be produced in perturbative gluon splittings, 
g — > bb, but this is estimated to be a very small contribution at LEP1 energies [[15] . 
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In ref it is assumed that the next qq pair in the string breakup has the same 
probability to be a uu as a dd pair. This is not necessarily the case, even if the fragmentation 
is iso-spin invariant. Consider a toy model in which only pions are produced. Then iso-spin 
invariance implies equal numbers of tt + , tt°, and tc~ , which means that a dd pair must be 
followed by a uu with probability 2/3 and a dd only with probability 1/3. (Equal probabilities 
would give twice as many tt° as tt + or 7r~.) This is the type of flavour correlations discussed 
in ref ||. 

From our estimates the effect of a realistic degree of flavour correlations is small. It 
could conceivably affect an estimate of Bn resonance by distorting the background. The 
flavour correlations imply that the excess positive charge is over compensated in the first 
rank meson, followed by a negative contribution in the second rank, a still smaller positive 
contribution in the third rank meson, etc. We want in the future to estimate the magnitude 
of such a possible distortion, to see if it could have a noticeable affect on the experimental 
results. We note that the OPAL Collaboration estimates that (Btc) resonances contributes 
with at least 20% of the total production of B mesons in Z° — > bb fragmentation |14 . 



In [0 a mass estimate for the (B7r) resonance is obtained from extrapolation of charmed 
meson data to B mesons. The estimate puts an upper limit on the resonance mass at 5.8GeV. 
This estimate is supported by OPAL data where an excess, as compared with simulations 
where no resonance meson production was included, of B + 7r pairs is found at 5.7GeV. The 
(B7r) resonance will eventually decay to a Btc pair close in phase space, and contribute to 
this excess. 

Pions produced in the center of the string have a fairly large momentum relative to 
the B° meson. If we want to retain B°7r pairs close in phase space, the relative momentum 
between central pions and the B puts kinematical constraints on the B energy, forcing the B 
energy to be small. We have seen above that the average xe for B mesons is 0.70, so central 
pions will have small chances of making a small mass pair with the B. Low xe mesons can 

The method to trace CP violating B mesons proposed in ref depends on the corre- 
lation of B° (B°) mesons with pions nearby in the phase space, and the detection of a CP 
eigenstate f as a B° or B° decay product in conjunction with a similar pion. 

The suggested method leads to a dilution of any real CP violating asymmetry. Starting 
with the time integrated asymmetry |16| 



, m _ r(Bj ^ f) - r(B° t=0 ^ f) 
Uj "r(Bto-f) + r(B° t=0 ^f)' [b) 

where f denotes a CP eigenstate and r(B° =0 — > f) is the decay width when the state f is 
produced decay product of a state starting out clS cl B° at t=0. 

Expressing the asymmetry in final states we get 

= 1 [N(Tn-) + N(Tn~)}N i+ - [iy(Tvr+) + jV(T7r+)]iV f - 
U ' 1 + x 2 [JV(Ttt-) - iV(T7r-)]A^ f+ + [N{Tn+) - N(Tn+)}N { - ' 1 ' 

T is a flavoured state we know to come from a B°, and is used to determine the flavour of the 
B meson. iV(T7r) is the relative number (probability) of states T in conjunction with a pion. 
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Figure 4: Invariant mass multiplicity distributions. 

a) Events where all possible B°7i pairs are considered. The full lines are data from events 
without any tensor mesons. The upper full line is for B°ti + pairs, and the lower full line is 
for B°ti- pairs. The dashed lines are for events with a tensor meson fraction of 30%. Again 
the upper (dashed) line is for B°tt + and the lower for B°n~ . 

b ) Events where only the fastest pion is paired with the B meson. The different lines describe 
the same sort of data as in a). 



Nft denotes the number of final states f7r . x is the mass mixing parameter x = (Am/T). 
(Notation taken from []J.) 

Considering a charge symmetric production process (pp or e + e~) where B° and B° 
production should be equal, ie. iV(T7r~) = iV(T7r + ) and iV(T7r + ) = iV(T7r~), gives 

Aobs ~ N {+ + N { - ~ iV(Tvr+) + iV(TTT-) { +X } U h {7) 

and we see that the tagging process suppresses the experimentally observed asymmetry. 
We are going to estimate the efficiency of the suggested tagging method by calculating the 
dilution factor, 

iV(T7r + ) -JV(Ttt-) , , 



N(Trr+) + N(Tir~)' 



of Eq (0) for the reaction e + e~— > bb — > X as a function of the B7r pair mass. From the B 
production considerations above we expect an excess of low mass B°7r + pairs giving a non 
vanishing dilution factor. 

In the calculations we use Jetset with the softer fragmentation distribution, Eq @, 
since this is in best agreement with experimental data. As the tagged state we use the B° 
meson and accept all possible B°7r pairs. In Fig f|a we show the invariant mass multiplic- 
ity distribution for B°7r pairs from 10 6 mixed flavour events at LEP1 energies. Of these 
approximately 218000 events will be bb events giving b hadrons. In the figure we clearly 
see a difference between the B°7r + and B°7r~ distributions. Continuing to large m(Bn) the 
difference decreases since the B7r-charge correlations vanishes. 
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Figure 5: Dilution factor as a function of invariant mass. 

a) Events where all possible B°n pairs are considered. The full line is for events without any 
tensor mesons, and the dashed line is for events with a tensor meson fraction of 30%. 
h) Events where only the fastest pion is paired with the B meson. The two lines describe the 
same sort of data as in a). 



In the calculation of the average dilution factor we only include Bit pairs with masses 
below 5.8GeV. This restriction eliminates the chance of using central pions in the pair and 
emphasize the contributions from (B71") resonances which are more probable to give a right 
charged pion. Accepting all possible pairs with m(B°7r) < 5.8GeV the average dilution factor 
is £ =0.13-0.17. The variation in £ is mainly due to different production fractions of tensor 
mesons allowed in the simulations, and we note that a larger fraction of resonance meson 
production will enhance the efficiency of the tagging method. 

The efficiency of the proposed tagging method will increase if only the most correlated 
pion in an event is taken into account. This pion should be the fastest one in the event. 
In Fig |2| we have plotted the rapidity difference between rank 2 and rank 3 mesons in the 
Lund string. Indeed we see that in general the fastest meson is the rank 2 particle, ie. the 
neighbour to the b hadron. (In Fig |^ we show distributions for p mesons, but similar results 
are obtained for any directly produced particles.) Considering only this pion in every event 
gives a more peaked invariant mass distribution as can be seen in Fig ^b. The average 
dilution factor, for pairs with m(B7r) < 5.8GeV, increases to 0.17-0.20. Again the dilution 
factor depends on the fraction of tensor mesons produced in the simulations. 

The distributions in Fig f| can be used to calculate the dependence of the dilution factor 
on the B7r invariant mass. The result is plotted in Fig We see that at m(B7r) = 5.5GeV 
the dilution factor is about 0.10 and grows to a maximum at m(B7r) = 5.8GeV. The growth 
is much larger if only the fastest pion is used as can be seen in Fig [5]b. The large difference 
in the maximum of £ = around 5.8GeV depend on whether tensor mesons are produced 
or not, since the decay of B** mesons contributes to the B°7r + and therefore enhance the 
dilution factor. In Fig |5|b £ remains around 0.25 above m(B7r) = 6GeV, but very low rates 
are expected above this mass, cf. Fig |]b. 

When the tensor meson production rate is assumed to be within 20-30%, we calculate 
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the dilution factor in the invariant mass range 5.6 < m (b°tt) < 5.9GeV to be 0.31-0.35. 
The OPAL Collaboration finds a similar value, 0.364±0.024, in their B + 7r studies |T3]], where 
it is expected that at least 20% of the B mesons are produced as decay products of heavier 
resonances at LEP1 

We have also performed the same calculations with an MC where we included flavour 
correlations between neighbouring qq pairs. The results from these runs are very similar to 
those found above without such flavour correlations. 



4 Conclusions 

A modification of the original Lund string fragmentation function (Eq (0)), where the re- 
duction in the "colour coherence area" due to the large mass of heavy quarks is taken into 
account, will give a softer energy distribution for heavy mesons. Experimental data are re- 
produced by tuning the parameter tq to 1.05±0.07 (cf. Eq (|3|)), which is consistent with the 
theoretically derived Tq = 1. The Lund string model, together with the modified splitting 
function, also gives rapidity distributions in agreement with PQCD assuming LPHD or clus- 
ter fragmentation. The Jetset MC is found to perform better when the modified splitting 
function is used, as compared to the the original splitting function or the Peterson et al. 
fragmentation model. 

The possibility to experimentally observe CP violating B mesons using pion charge to 
tag the flavour of the leading B° is always diluted. This dilution depends on the fraction of 
resonance B meson production. In the OPAL measurement fbl it is found that at least 20% 



of the B mesons are produced as decay products of heavier resonances at LEP1. We estimate 
the dilution factor, defined in Eq (||), in the invariant mass range, 5.6 < m (P ? 1 ") < 5.9GeV 
to be 0.31-0.35. The uncertainty is due to the fraction of resonance production, here we 
have used 20-30%. The OPAL Collaboration finds a similar value, 0.364±0.024, in their 



B 7r studies [14 
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